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ABSTRACT 

We report on Chandra observations of a sample of 11 optically luminous (Mb< —28.5) quasars at 
z=3.96^.55 selected from the Palomar Digital Sky Survey and the Automatic Plate Measuring Facility Survey. 
These are among the most luminous z ^ 4 quasars known and hence represent ideal witnesses of the end of 
the “dark age ”. Nine quasars are detected by Chandra, with « 2-57 counts in the observed 0.5-8 keV band. 
These detections increase the number of X-ray detected AGN at z ^ 4 to ~ 90; overall, Chandra has detected 
« 85% of the high-redshift quasars observed with snapshot (few kilosecond) observations. PSS 1506-1-5220, 
one of the two X-ray undetected quasars, displays a number of notable features in its rest-frame ultraviolet 
spectrum, the most prominent being broad, deep Si IV and C IV absorption lines. The average optical-to-X-ray 
spectral index for the present sample ((aox)=—L88±0.05) is steeper than that typically found for z ^ 4 quasars 
but consistent with the expected value from the known dependence of this spectral index on quasar luminosity. 

We present joint X-ray spectral fitting for a sample of 48 radio-quiet quasars in the redshift range 3.99-6.28 
for which Chandra observations are available. The X-ray spectrum (« 870 counts) is well parameterized by a 
power law with r=L93^QQg in the rest-frame « 2^0 keV band, and a tight upper limit of Ah sa5xl0^' cm ^ 
is obtained on any average intrinsic X-ray absorption. There is no indication of any significant evolution in the 
X-ray properties of quasars between redshifts zero and six, suggesting that the physical processes of accretion 
onto massive black holes have not changed over the bulk of cosmic time. 

Subject headings: galaxies: active — galaxies: nuclei — quasars: general — X-rays: galaxies 


1. INTRODUCTION 

Over the last few years. X-ray observations have pro¬ 
vided important information about the population of luminous 
quasars populating the Universe at z « 4-6, i.e., at the end of 
the “dark age” (e.g., Rees 1999). Since the work by Kaspi, 
Brandt, & Schneider (2000) with archival ROSAT observa¬ 
tions, X-ray studies of high-redshift quasars have taken ad¬ 
vantage of the capabilities of the new generation of X-ray 
observatories: Chandra and XMM-Newton. In particular, 
Chandra snapshot observations of z ^ 4 Active Galactic Nu¬ 
clei (AGN) have proven to be effective in providing basic 
X-ray information, such as X-ray fluxes, luminosities, and 
broad-band properties, on some of the most luminous ob¬ 
jects in the Universe (Vignali et al. 2001a, 2002, 2003a,b, 
hereafter VOla, V02, V03a, and V03b, respectively; Brandt 
et al. 2002; Bechtold et al. 2003; Bassett et al. 2004). 
Using these observations, it has also been possible to col¬ 
lect enough source counts to allow joint X-ray spectral fit¬ 
ting (V03a; V03b; Vignali, Brandt, & Schneider 2004). This 
kind of X-ray study is enabled by the sharp on-axis PSF and 
the low instrumental background of Chandra. On the other 
hand XMM-Newton, with its larger effective area, has al¬ 
lowed individual X-ray spectral htting for a handful of the 
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most X-ray luminous z > 4 quasars with medium-to-long ex¬ 
posures (Ferrero & Brinkmann 2003; Grupe et al. 2004; Far- 
rah et al. 2004; see also the Chandra spectral results of the 
z=5.99 SDSS 130608.26-1-035626.3 quasar reported recently 
by Schwartz & Virani 2004). 

It has therefore become possible to provide an overall pic¬ 
ture of the radio-quiet and radio-loud quasar populations at 
the highest redshifts. However, X-ray studies of the most op¬ 
tically luminous quasars at z > 4 have been conducted for a 
limited number of objects, the majority of which have been se¬ 
lected from the Digital Palomar Sky Survey (PSS; e.g., Djor- 
govski et al. 1998) and presented by V03a. The goal of the 
present study is to extend previous X-ray investigations of the 
properties of luminous quasars at z ^ 4 using Chandra snap¬ 
shot observations; the target sample comprises 11 quasars in 
the redshift range z=3.96^.55 (Mb~ —28.5 to —29.5) se¬ 
lected from the PSS and the Automatic Plate Measuring Fa¬ 
cility Survey (BRI; e.g., Irwin, McMahon, & Hazard 1991; 
Storrie-Lombardi et al. 1996,2001). 

Throughout the paper we adopt Hq=10 km s^' Mpc^*, 
f2M=0.3, and fiA=0.7 (e.g., Spergel et al. 2003). 

2. Chandra OBSERVATIONS AND DATA REDUCTION 

All of the quasars presented here were targeted by Chandra 
during Cycle 4 with snapshot (ft! 4.0-5.1 ks) observations. 
The observation log is shown in Table 1. All of the sources 
were observed with the Advanced CCD Imaging Spectrom¬ 
eter (ACIS; Garmire et al. 2003) with the S3 CCD at the 
aimpoint. Faint mode was used for the event telemetry for¬ 
mat, and ASCA grade 0, 2, 3, 4 and 6 events were used in 
the analysis, which has been carried out using the Chandra 
Interactive Analysis of Observations (CIAO) Version 2.3 soft¬ 
ware. No background flares are present in these observations. 
Source detection was carried out with WAVDETECT (Freeman 
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et al. 2002) similarly to VOla, V03a, and V03b, using wavelet 
transforms (with wavelet scale sizes of 1, 1.4, 2, 2.8, and 
4 pixels) and a false-positive probability threshold of 10^^. 
Given the small number of pixels being searched due to the 
known source positions and the sub-arcsec on-axis angular 
resolution of Chandra, the probability of spurious detections 
is extremely low; most of the sources were in fact detected at 
a false-positive probability threshold of 10^®. 

The X-ray counts detected in the ultra-soft band 
(0.3-0.5 keV), the soft band (0.5-2 keV), the hard band 
(2-8 keV), and the full band (0.5-8 keV) are reported in Ta¬ 
ble 2. The counts have been derived using WAVDETECT and 
checked with manual aperture photometry. Nine quasars have 
been detected, with full-band counts ranging between 2 and 
57. The X-ray positions of the detected sources are within 
0.1-0.7" of the optical positions; this is consistent with the 
expected Chandra ACIS positional error. 

In Table 2 the band ratios and the effective photon indices 
for the X-ray detected sources are also reported. Both of these 
quantities have been corrected for the quantum efficiency de¬ 
cay of ACIS at low energies, caused by molecular contami¬ 
nation of the ACIS filters, using a time-dependent correction. 
To apply the correction, we adopted the ACISABS model (see, 
e.g., Chartas et al. 2002).® We note that, given the dates of 
our observations, the differences between the adoption of ei¬ 
ther this model or that proposed by the Chandra X-ray Cen¬ 
ter (CXC),^ available with the latest CIAO release, are small; 
therefore, for consistency with our previous work (e.g., V03a; 
V03b), we will adopt the ACISABS model throughout this pa¬ 
per. 

The full-band images of eight detected quasars, adaptively 
smoothed using the algorithm of Ebeling, White, & Rangara- 
jan (2005), are shown in Eigure 1. Eor BR 1117—1329, 
a known broad absorption-line quasar (BALQSO; Storrie- 
Lombardi et al. 1996) from which two counts have been de¬ 
tected (see §2.1 for discussion), and the two non-detections, 
the raw images are shown. Overall, in the regions of pro¬ 
jected distance « 650 x 650 kpc^ centered on the quasar po¬ 
sitions, we did not find any significant excess of companions 
with respect to the cumulative number counts from X-ray sur¬ 
veys (e.g., Bauer et al. 2004); this is in agreement with previ¬ 
ous X-ray searches around high-redshift quasars (e.g., V03a; 
V03b). Similarly, the absence of X-ray extension for all of 
our target quasars suggests that significant gravitational lens- 
ing (e.g., Comerford, Haiman, & Schaye 2002) is not present. 

Although the observations are short, we have searched for 
X-ray variability by analyzing the photon arrival times in the 
full band using the Kolmogorov-Smirnov test. No significant 
variability was detected from any of our quasars. This is not 
surprising given the small numbers of detected counts and the 
w 15-minute rest-frame exposure times. 

2.1. Two-Count Sources: BR 1117—1329 and 
PSS 1506+5220 

As was the case for the z=5.21 quasar SDSS 1208-1-0010 
(VOla) and the z=4.80 quasar SDSS 0756-1-4502 (V03b), 
BR 1117—1329 can also be considered a secure detec¬ 
tion with only two counts (in the same pixel) in the 
0.5-2 keV band. In fact, the source is detected by WAVDE¬ 
TECT using a false-positive threshold of 10^^. Monte-Carlo 

^ See http://www.astro.psu.edu/users/chartas/xcontdir/xcont.html 

^ For further details, see http://asc.harvard.edu/ciao/why/acisqedeg.html 


simulations (see §2.1.1 of V03b for details) indicate that the 
significance of this detection is 4c in the soft band. 

In contrast, PSS 1506-1-5220 cannot be considered detected 
(according to a Monte-Carlo analysis), although two full-band 
counts are present within the 2"-radius circle centered on its 
optical position (see Eig. 1). 

3. HOBBY-EBERLY TELESCOPE OBSERVATIONS 

We conducted imaging and spectroscopic observations of 
several of the quasars using the 9-m Hobby-Eberly Telescope 
(HET; Ramsey et al. 1998). These data were taken within 
1-10 months of the X-ray observations (< 8 weeks in the 
quasar rest frames) to minimize the effects of variability in 
the X-ray/optical comparisons. Variability can introduce sig¬ 
nificant uncertainties into estimates of key parameters (e.g., 
the optical-to-X-ray spectral index; see §4). The observations 
were all obtained with the Marcario Low-Resolution Spectro¬ 
graph (LRS; Hill et al 1998a,b; Cobos Duenas et al. 1998). 

3.1. Photometric Observations 

Two-minute R-band LRS images were acquired for four 
of the 11 sources (PSS 0747-f4434, PSS 15064-5220, 
PSS 16464-5514, and PSS 23444-0342) during the fall-winter 
2003^. Using the k,A' xA’ LRS images and finding charts 
of the fields, it was possible to obtain approximate photomet¬ 
ric calibrations using the stellar objects in the fields. Although 
the level of uncertainty in these comparisons is substantial, we 
did not detect any significant long-term variations (i.e., more 
than 50%) in the quasar brightnesses from their published val¬ 
ues. To have a uniform analysis, we have adopted the APM 
magnitudes^ to derive the optical and broad-band parameters 
of the quasars. 

3.2. Optical Spectrum of PSS 1506+5220 

Two of the quasars, PSS 1058-41245 and PSS 1506-45220, 
do not have published spectra. Immediately after acquiring 
the R-band image of PSS 1506-45220 on January 19, 2004 
(described in §3.1), we obtained a 15-minute spectrum of the 
object using the LRS. The spectrum, which is displayed in 
Eig. 2, was acquired using a 2" slit, an OG515 blocking filter, 
and a 600 line mm^^ grism; this configuration produced a 
spectrum covering the wavelength range 6300-9100 A at a 
resolving power of « 1100. 

The spectrum displays a number of notable and complex 
features, the most prominent being broad, deep Si IV and CIV 
absorption lines between redshifts of about 3.68 and 3.84 (a 
velocity width of approximately 10,000 km s^^). The sharp 
rise at the blue end of the spectrum is the red wing of the Ly¬ 
man a emission line. Also present is a strong Mg II absorption 
doublet (total rest equivalent width of nearly 7 A). 

The previously reported redshift of PSS 1506-45220 is 
4.18;^ however, accurate redshifts of z > 4 BALQSOs are dif¬ 
ficult to determine. In Pig. 2 we have indicated the expected 
position of the C IV emission line; no feature is found at this 
location. There is a suggestion of a broad emission feature 
centered at 7870 A; if this is the C IV line, the redshift would 
be approximately 4.08 (there is a hint of a narrow absorption 
feature in the midst of this putative emission line). 

G. Djorgovski has kindly showed an unpublished Keck 
spectrum of this object to us; his data extend far below the Ly¬ 
man a emission line. If the peak of this line is associated with 

* See http://www.ast.cam.ac.uk/~mike/apmcat/interface.html 

^ See http://www.astro.caltech.edu/~george/z4.qsos 
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the rest wavelength of Lyman a, then the redshift is indeed 
4.18. Given the significant distortions in the lines produced 
by the absorption troughs, we conclude that it is presently im¬ 
possible to assign a redshift with an accuracy of better than 
w 0.1 for this object. The ejection velocities of the BAL fea¬ 
tures would range between 20,000 km s^ * and 25,000 km s^ *, 
depending on the quasar redshift. Throughout the analysis in 
this paper, we will adopt the published value of 4.18 for the 
redshift of this object (see Footnote 9). 

The strong Mg II absorption feature suggests the presence 
of a massive galaxy at a redshift of 1.4711 along the line-of- 
sight to the quasar. Given this result, the high luminosity of 
PSS 1506-1-5220 may partially arise from amplification via 
gravitational lensing. 

4. MULTI-WAVELENGTH PROPERTIES OF THE SAMPLE 

The principal optical. X-ray, and radio properties of our 
target quasars are listed in Table 3: 

Column (1). — The name of the source. 

Column (2). — The Galactic column density (from Dickey & 
Lockman 1990) in units of lO^*' cm^^. 

Column (3). — The monochromatic rest-frame AB 1450 ( 1 + 2 ) 
magnitude (ABi450(i+2)=-2.5 log /1450 a - 48.6; Oke 
& Gunn 1983). The magnitudes reported in Table 3 
have been derived from the ARM B-band magnitudes 
(corrected for Galactic extinction; Schlegel, Finkbeiner, 
& Davis 1998) assuming the empirical relationship 
ABi45o{i+z)=^“0.684z-|- 3.10, which is effective in the 
redshift range covered by our sample. 

Columns (4) and (5). — The 2500 A rest-frame flux 

density and luminosity. These were computed from the 
AB 1450 ( 1 + 2 ) magnitude assuming an optical power-law slope 
of a = —0.79 (Bv v“; Fan et al. 2001) to allow direct 
comparison with the results presented in VOla, V03a, and 
V03b. Note that the 2500 A rest-frame flux densities and 
luminosities are reduced by 15% when the power-law 
slope of the optical continuum is changed to a = —0.5 (e.g., 
Schneider et al. 2001; Vanden Berk et al. 2001). 

Column (6). — The absolute B-band magnitude computed 
using a = —0.79. When a = —0.5 is adopted for the extrap¬ 
olation, the absolute B-band magnitudes are fainter by « 0.35 
mag. 

Columns (7) and (8). — The observed count rate in the 
0.5-2 keV band and the corresponding flux, corrected for 
Galactic absorption and the quantum efficiency decay of 
Chandra ACIS at low energy. The fluxes have been calculated 
using PIMMS (Mukai 2002) and a power-law model with 
r=2.0; this is a good parameterization of samples of z « 0-3 
quasars (e.g., George et al. 2000; Reeves & Turner 2000; 
Page et al. 2003; Piconcelli et al. 2005) and z > 4 quasars 
(V03a; V03b; see § 6 ). 

Columns (9) and (10). — The rest-frame 2 keV flux density 
and luminosity, computed assuming r=2.0 and corrected 
for the quantum efficiency decay of Chandra ACIS at low 
energy. 

Column (11). — The 2-10 keV rest-frame luminosity. 
Column (12). — The optical-to-X-ray power-law slope, ttox, 
defined as 

_ log(/2kev//2500 a) /ix 

“ log(V2kev/V2500A) 

where /2 keV and /2500 a are the rest-frame flux densities at 
2 keV and 2500 A, respectively. The « la errors on ttox have 
been computed following the “numerical method” described 


in § 1.7.3 of Lyons (1991) considering both the statistical un¬ 
certainties on the X-ray count rates and the effects of reason¬ 
able changes in the X-ray and optical continuum shapes (see 
§ 3 of VO 1 a for details). 

Column (13). — The radio-loudness parameter (e.g., Keller- 
mann et al. 1989), defined as R=fs ghz// 44 oo A frame). 
The rest-frame 5 GHz flux density was computed from the 
FIRST (Becker, White, & Helfand 1995), NVSS (Condon 
et al. 1998) or Carilli et al. (2001) observed 1.4 GHz flux 
density assuming a radio power-law slope of a = —0.8. The 
upper limits are at the 3a level. The rest-frame 4400 A flux 
density was computed from the ABi 45 o(i+ 2 ) magnitude assum¬ 
ing an optical power-law slope of a = —0.79. Typical radio¬ 
loudness values are > 100 for radio-loud quasars (RLQs) and 
< 10 for radio-quiet quasars (RQQs). 

For the two quasars with the lowest declinations, 
BR 0418—5723 and BR 2213—6729, no sensitive radio data 
are available. None of the remaining quasars is radio loud. 

5. X-RAY PROPERTIES OF THE SAMPLE 

In our previous study of a sample of optically luminous 
(Mb« -28.4 to -30.2) z > 4 PSS quasars (V03a), we de¬ 
rived the basic X-ray properties for some of the most lumi¬ 
nous objects in the Universe and performed joint X-ray spec¬ 
tral fitting. No spectral evolution in quasar X-ray continua at 
high redshift was found. With the current sample of 11 PSS 
and BRI quasars, our goal is to improve the results presented 
in V03a. The new Chandra observations, coupled with some 
other observations retrieved from the archive,**' 

Most of the quasars presented in this paper (filled triangles 
in Fig. 3 for the nine X-ray detections and large down ward¬ 
pointing arrows for the two upper limits) populate the re¬ 
gion close to the aox=—1-80 slanted line (computed assuming 
z=4.3, the average redshift of the current sample). They have 
an average (aox)=— 1-88±0.05 (computed using the ASURV 
software package Rev 1.2; LaValley, Isobe, & Feigelson 
1992; the quoted errors represent the standard deviation of 
the mean) which is steeper than the average value found using 
the whole population of RQQs at high redshift with available 
X-ray information ((aox)=— U76±0.02). 

This finding is consistent with the fact that 
the current Chandra sample is highly luminous 
((log/25 oo)= 32.1 erg s^* Hz^*), and from previous work 
(e.g., Avni, Worrall, & Morgan 1995; Vignali, Brandt, & 
Schneider 2003, hereafter VBS03; see also Strateva et al., 
in preparation) ttox has been found to depend upon optical 
luminosity, i.e., more optically luminous quasars have steeper 
(more negative) ttox- The expected ttox for the current sample 
would be « —1.8 ±0.1. We note that ttox for BR 0353—3820 
is flat ( 0 Cox=~l- 54 /gQg), much below the average value for 
the whole sample; it is possible that this flat ttox is due 

See http://www.astro.psu.edu/users/nieI/papers/highz-xray-detected.dat 
for a regularly updated compilation of X-ray detections and tight upper limits 
at 4 > 4. Details about the analyzed Chandra archival observations are re¬ 
ported in the Appendix, have increased the number of X-ray detected AGN 
at 4 > 4 (mostly quasars) to « 90. Overall, Chandra has detected « 85% 
of the high-redshift quasars observed with snapshot observations. The pro¬ 
gram of snapshot observations has allowed filling of the soft X-ray flux versus 
A 6 i 45 o(i+z) plot shown in Eig. 3 with many X-ray detections over the last few 
years. This plot is effective in showing the large advances made recently in 
the knowledge of the properties of high-redshift quasars (only RQQs are re¬ 
ported); for compaiison, note the paucity of data points in Fig. 2 of Brandt 
et al. (2001). Most of the upper limits shown in Fig. 3 have been obtained 
using archival ROSAT observations; some of these are reported here for the 
first time (see Appendix and Table Al for details). 
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to some source variability occurring over the time interval 
between the optical and X-ray observations. 

As discussed in previous works (e.g., VOla; V03b), there 
is a correlation between ABi 45 o(i+j) and 0.5-2 keV flux (see 
Fig. 3). To have a homogeneous sample, we have used 
only the z ^ 4 optically selected RQQs observed by Chan¬ 
dra', the significance of the correlation has been evaluated 
using the methods within ASURV, in particular the general¬ 
ized Kendall’s x (Brown, Hollander, & Korwar 1974). The 
ABi 45 o(i_|_j) vs. 0.5-2 keV flux correlation is significant at the 
> 99.9% level; according to the Estimate and Maximize (EM; 
Dempster, Laird, & Rubin 1977) regression algorithm, it can 
be parameterized by 

log (Ex/erg cm-2 s^') = [-(0.293 ±0.057) 
xABi450(i+A-(8-786 ±1.059)], 

plotted as a dotted line in Fig. 3. 

Figure 3 also highlights the large discrepancy in our knowl¬ 
edge of the properties of the high-redshift AGN population 
in its entirety. While X-ray information is presently avail¬ 
able for large numbers of optically luminous quasars at z ^ 4 
(i.e., those located on the right-hand side of Fig. 3 and plotted 
as filled symbols and downward-pointing arrows), only a few 
AGN at fainter optical magnitudes are detected in the X-rays 
(open squares). These have all been discovered by moderately 
deep and deep X-ray surveys with ROSAT (Schneider et al. 
1998) and Chandra (Silverman et al. 2002; V02; Castander 
et al. 2003; Treister et al. 2004). 

Two quasars in the present sample have not been de¬ 
tected by Chandra: PSS 1506±5220 (discussed in §2.1) 
and PSS 2344±0342 (no X-ray counts within a 2"-radius 
circle). There are two possible explanations for the non¬ 
detections of these two quasars in the X-ray band; they are 
either strongly absorbed or intrinsically X-ray weak. If their 
non-detection is due to the presence of intrinsic X-ray ab¬ 
sorption, then column densities of Ah ± 1.5 x 10^^ cm^^ and 
Ah 3.3 X 10^^ cm^^, respectively, are required to repro¬ 
duce the Chandra constraints (using the optical-to-X-ray cor¬ 
relation of VBS03 and assuming a photon index r=2.0). It is 
interesting to note that the HET spectrum of PSS 1506±5220 
displays a number of notable and complex features, the most 
prominent being broad, deep Si IV and C IV absorption lines. 
For the other X-ray undetected quasar of the present sample, 
PSS 2344±0342, the optical spectrum shows two very high 
column density damped Lyman a absorption (DLA) systems 
(Peroux et al. 2001; Prochaska et al. 2003); however, these 
are unlikely to be responsible for all of the absorption sug¬ 
gested by the X-ray constraints (see, e.g., Elvis et al. 1994 for 
discussion on the expected X-ray absorption by DLAs). 

We also note that the X-ray weak BR 1117—1329 (detected 
with two counts, see §2.1 and Table 2) shows blueshifted 
broad absorption lines for O VI, N V, Si IV, and C IV (Storrie- 
Lombardi et al. 1996). 

6 . JOINT SPECTRAL FITTING: A QUASAR X-RAY SPECTRAL 
TEMPLATE AT Z « 4-6.3 

Previous studies have provided basic information on the 
X-ray spectral properties of luminous, optically selected z > 4 
RQQs, in particular on their photon indices and column den¬ 
sities. These results have been achieved either with direct 
X-ray spectral fitting of a few X-ray luminous quasar spec¬ 
tra observed with XMM-Aewfon (e.g., Ferrero & Brinkmann 
2003; Grupe et al. 2004; Farrah et al. 2004) and Chandra 


(Schwartz & Virani 2004) or using joint X-ray spectral fit¬ 
ting analysis, i.e., simultaneous fitting of sizable numbers of 
Chandra spectra of z > 4 quasars (e.g., V03a; V03b; Vig- 
nali et al. 2004). Since these objects are among the most 
luminous quasars at z > 4, they do not represent the major¬ 
ity of the high-redshift quasar population. Similarly to low- 
and intermediate-redshift RQQs (e.g., George et al. 2000; 
Reeves & Turner 2000; Page et al. 2003; Piconcelli et al. 
2005), high-redshift RQQs are, on average, characterized by 
power-law continua with F « 1.8-2.0, with no indication for 
widespread absorption. Joint X-ray spectral fitting studies 
are important given the difficulty in obtaining quality X-ray 
spectra of high-redshift quasars using reasonable exposures 
(« 40-60 ks) with XMM-Newton. Furthermore, combining 
the spectra of many high-redshift quasars is a practical way to 
derive their average X-ray spectral properties; in fact, studies 
limited to a few X-ray luminous quasars may not be represen¬ 
tative of the high-redshift quasar population as a whole. 

To provide further constraints on the average X-ray proper¬ 
ties of high-redshift quasars, we have performed joint fitting 
using all the available Chandra data for RQQs in the redshift 
range z=3.99-6.28. Neither X-ray selected quasars (see V02 
for X-ray spectral analyses of X-ray selected z > 4 AGN in the 
Chandra Deep Field-North) nor radio-selected quasars (with 
the exception of FIRST 0747±2739, which is radio quiet; see 
Bassett et al. 2004) have been considered in the following 
analysis. Source counts have been extracted for each quasar 
using a circular region of 2" radius centered on the X-ray cen¬ 
troid of the quasar. Although background is typically negli¬ 
gible in these Chandra observations, we have extracted back¬ 
ground counts using annuli of different sizes (to avoid con¬ 
tamination from nearby X-ray sources) centered on the quasar 
position. Only sources with > 2 full-band counts have been 
considered in the analysis. The selection criteria adopted in 
this study (mainly the RQQ selection and the fact that only 
Chandra data are used) reduced the number of z ^ 4 quasars 
with usable X-ray information from « 90 to 48, for a total ex¬ 
posure time of 244.1 ks. We have w 870 net source counts in 
the observed 0.3-8 keV band, corresponding to the rest-frame 
« 1.5-58 keV band; the average number of counts per object 
is « 18. Figure 4 summarizes the principal properties of the 
sample used in this X-ray spectral study. The average redshift 
of the sample is z=4.57, while its median is z=4.43. 

Spectral analysis was carried out with XSPEC Version 
11.3.0 (Arnaud 1996) using unbinned data and the C-statistic 
(Cash 1979); this statistical approach allows one to retain all 
spectral information and associate with each quasar its own 
Galactic absorption column density and redshift (for the fit¬ 
ting with intrinsic absorption, see below). Errors are quoted 
at the 90% confidence level for one interesting parameter 
(AC = 2.71; Avni 1976; Cash 1979), unless stated otherwise. 
Solar abundances (from Anders & Grevesse 1989) have been 
adopted for the absorbing material, although optical studies 
indicate that high-redshift quasars are often characterized by 
supersolar abundances of heavy elements (e.g., Hamann & 
Ferland 1999; Dietrich et al. 2003a,b). Note that doubling 
the abundances in the fit produces a reduction in the column 
density of a factor of w 2. 

We started the joint X-ray spectral fitting using a power- 
law model, leaving the power-law normalizations for all of 
the quasars free to vary (to account for the different fluxes 
of the sources), plus Galactic absorption. The lack of sig¬ 
nificant data-to-model residuals for the fit suggests that the 
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assumed model is acceptable; the resulting photon index, 
r=1.93±0.09, is consistent with previous estimates for z > 4 
quasars. The addition of neutral intrinsic absorption at each 
source’s redshift does not improve the fit significantly; there 
is no evidence for widespread absorption toward the lines- 
of-sight of the quasars under investigation, the column den¬ 
sity (i.e., the counts-weighted average column density) being 

5.0 X 10^* cm^^ (see the solid contours in Fig. 5a showing 
the 68, 90, and 99% confidence levels for the photon index vs. 
column density; see also Table 4). If our targets had intrinsic 
absorption of Nh ~ (2-5) x 10^^ cm^^, as found in two RQQs 
at z « 2 using ASCA data (Reeves & Turner 2000), we would 
be able to detect it given the upper limit on the column density 
we derived. The spectral results do not change significantly if 
only the X-ray data in the rest-frame interval common to all of 
the quasars of the sample (w 2.2^0 keV band; « 840 counts) 
are taken into account (see the dashed contours in Fig. 5a). In 
this case, the photon index is 1.95 ±0.10, while the upper 
limit on the column density remains essentially unchanged. 

The sample does not appear to be biased by the presence 
of a few high signal-to-noise ratio objects." This has been 
checked during the X-ray spectral analysis by removing some 
of the quasars with higher photon counting statistics. In par¬ 
ticular, the removal of the three sources with ^60 counts 
(BR 0353-3820, PSS 0926±3055, and PSS 1326±0743) 
leaves the photon index unchanged (r=1.92±0.10), but the 
constraint on the column density intrinsic to the quasars is less 
tight ( 7.5 X 10^* cm see Table 4). We also performed 

joint spectral fitting using the most and least optically lumi¬ 
nous RQQs of the sample as well as the highest and lowest 
redshift RQQs; the results of the fitting are reported in Table 4. 
Within the statistical uncertainties due to the limited number 
of X-ray counts, no evidence for a significant dependence of 
the X-ray spectral parameters upon either the luminosity or 
the redshift of the quasar subsample has been found. 

Given the limited counting statistics available and the fact 
that a power-law model seems to be acceptable, we did not 
investigate more complex spectral models. The typical up¬ 
per limits on the equivalent width (EW) of a narrow, either 
neutral or ionized, iron Ka line are not stringent, lying in the 
range « 270-370 eV. Using the parameterization of the anti¬ 
correlation between the EW of the iron Ka line and source 
X-ray luminosity (the “X-ray Baldwin effect”; see Iwasawa 
& Taniguchi 1993 and Nandra et al. 1997) recently reported 
by Page et al. (2004), we would expect EW« 50 eV for our 
sample of RQQs (their average rest-frame 2-10 keV luminos¬ 
ity is 10^^ ^ erg s^*). The expected EW is below the upper 
limits provided by joint X-ray spectral fitting. 

To investigate the reliability of the joint X-ray spectral fit¬ 
ting results, we performed two consistency checks. First, we 
compared the photon index obtained from the joint spectral 
fitting with that derived from band-ratio analysis (i.e., the ratio 
of the 2-8 keV to 0.5-2 keV counts, corrected for the ACIS 
quantum efficiency decay); for the Galactic column density, 
we assumed the average value for the 48 quasars weighted 
by the exposure time of each observation. The resulting pho¬ 
ton index, r=1.89±0.08 (the uncertainty quoted here corre¬ 
sponds to the lo level, see Gehrels 1986 and the caption of 
Table 2), is consistent with our previous value. Secondly, 
we added all of the source and background counts using the 

*' For this reason, we have considered the ai 8.1 ks Chandra observation 
of SDSS 130608.26+035626.3 (Brandt et al. 2002) instead of the longer 
(Ri 120 ks) observation presented recently by Schwartz & Virani (2004). 


FTOOL task MATHPHA; similarly, redistribution matrix files 
(RMFs) and auxiliary response files (ARFs) were combined 
using ADDRMF and ADDARF, respectively, weighted by the 
individual exposure times. In this case the source counts 
were grouped into a spectrum such that each spectral bin con¬ 
tained at least 15 counts to allow fitting. Our fitting ob¬ 
tained r=1.97±0.11 and no evidence for intrinsic absorption 
(A^h .^S4.4 X 10^* cm ^), as for the previously adopted fit¬ 
ting procedures. The power-law model and the data-to-model 
residuals are shown in Fig. 5b. 

Fig. 6a shows the hard photon index (to avoid spectral com¬ 
plexities such as soft excesses at low energies) as a function of 
redshift for a compilation of optically selected RQQs. In par¬ 
ticular, at low (z < 1.8) redshift we have used the X-ray spec¬ 
tral information for the Bright Quasar Survey (PG; Schmidt & 
Green 1983) objects presented by Piconcelli et al. (2005; cir¬ 
cles); a few sources with Seyfert-like luminosity (Mb> —23) 
have been excluded from this analysis. At z « 1.8— 2.5, the 
RQQs from the Vignali et al. (2001) sample have been se¬ 
lected (diamonds). At the highest redshifts, the X-ray spectral 
results obtained over the last three years with either joint spec¬ 
tral fitting (triangles) or individual observations (squares) are 
shown. The low-redshift (z 1) quasars plotted in Fig. 6a are 
usually « 4 magnitudes less luminous than the z ^ 4 quasars 
under study; we are able to get down only to z ~ 1 with com¬ 
parably luminous objects. Although significant intrinsic scat¬ 
ter in the photon indices at all redshifts is present, there is no 
detectable systematic change in the photon-index distribution 
at high redshift. The apparently smaller dispersion in the pho¬ 
ton indices at high redshift than at lower redshift is at least 
partially due to the joint spectral fitting (used for some data 
points at z > 4) which tends to average any dispersion out. 
Although it is possible in principle that a redshift-dependence 
of r could be canceled out by a luminosity-dependence of T 
(which is not evident from Fig. 6b, where F is plotted against 
Mb), this would require somewhat of a conspiracy and cuts 
agains Occam’s razor. At present there is no overwhelming 
evidence that the photon index is significantly dependent upon 
quasar X-ray luminosity; for further details, see the contrast¬ 
ing results presented by Dai et al. (2004) and Reeves & Turner 
(2000). A partial correlation analysis of the photon index vs. 
redshift/X-ray luminosity for a large, uniformly selected sam¬ 
ple of quasars is required to address this issue properly. 

All of these spectral results confirm the lack of signifi¬ 
cant X-ray spectral evolution of luminous, optically selected 
quasars over cosmic time. Thus, despite the strong changes 
in environment and quasar number density that have occurred 
from z w 0-6, individual quasar X-ray emission regions seem 
to evolve relatively little (see Brandt et al. 2005). 

7. X-RAY UNDETECTED QUASARS: PECULIAR OBJECTS AT HIGH 
REDSHIET? 

Although the X-ray spectral results presented above are 
based on a significant fraction of the RQQs with Chandra 
observations, the selection of the sample might constitute a 
cause for concern, since we have excluded the X-ray faintest 
(with less than two counts) sources and non-detections. Over¬ 
all, these X-ray weak or undetected quasars represent « 25% 
of the RQQ population thus far observed by Chandra at z ^ 4. 
It is possible that in some cases the lack of an X-ray detec¬ 
tion is due to absorption, which should be ^ 10^^ cm^^ to 
reproduce the X-ray constraints (under reasonable assump¬ 
tions about the spectrum). This is an acceptable explanation 
for BALQSOs (40% of the the X-ray undetected quasars, see 
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Fig. 7). Past works have shown that BALQSOs are typically 
absorbed in the X-ray band (e.g., Brandt et al. 2000; Green et 
al. 2001; VOla; Gallagher et al. 2002). 

In the ttox histogram shown in Fig. 7, BAFQSOs pop¬ 
ulate preferentially the fainter part (i.e., more negative 
ttoxvalues) of the distribution. Although computed us¬ 
ing a limited number of objects, their average ttox (dot¬ 
ted line in Fig. 7) is —1.96 ±0.06 (—1.92 ±0.06 includ¬ 
ing SDSS 104433.04-012502.2 detected by XMM-Aewfon; 
see Brandt et al. 2001), while the optically selected non- 
BAF RQQs observed with Chandra have — 1.76 ± 0.02 (long- 
dashed line in Fig. 7; the median value is —1.77). It is also 
plausible, though not necessarily required by the present data, 
that some of the remaining X-ray undetected quasars at high 
redshift are characterized by absorption features which are not 
recognized because of either low spectral resolution (see, e.g., 
the case of SDSS 173744.87±582829.5 in V03b, where nar¬ 
row absorption features were found with the HEX) or limited 
spectral coverage (e.g., Maiolino et al. 2001, 2004; Goodrich 
et al. 2001). 

Given this, we believe that the criteria adopted in the selec¬ 
tion of the high-redshift non-BAF RQQs for X-ray spectral 
analysis (see §6) are not a cause for concern and the results 
we derived can be considered representative of the RQQ pop¬ 
ulation at z ^ 4. 

8. SUMMARY 

We have analyzed Chandra observations of a sample 
of 11 quasars at 3.96^.55 selected from the PSS and 
APM surveys. These are among the most luminous z ^ 4 
quasars known and represent ideal probes of the dawn of the 
modern Universe. Given the optical luminosities of these 
quasars (Mb~ —28.5 to —29.5), the present study can be 
considered complementary to past work (e.g., V03a; V03b). 
The principal results of this work are 

1. Nine of the quasars have been detected, thus increasing 
the number of X-ray detected AGN at z ^ 4 to w 90. In¬ 
cluding also the two X-ray non-detections, the average 
ttox of the present sample, — 1.88 ± 0.05, is steeper that 
typically found for z ^ 4 quasars, but consistent with 
the value expected applying the known dependence of 
this spectral index upon quasar luminosity. 

2. Two quasars are not detected, PSS 1506±5220 and 
PSS 2344±0342. If absorption intrinsic to the 
source were the cause of the X-ray weakness of 
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these two objects, then column densities larger than 
Ah ~ 1.5 X 10^^ cm^^ would be implied. In the case 
of PSS 1506±5220, we found evidence for notable 
features, the most prominent being broad absorption 
lines from Si IV and C IV. For this source it is possible 
that the UV absorber is linked to the absorption likely 
present in the X-ray band. This has been observed be¬ 
fore in other BAFQSOs at high redshift. 

3. Using all the available RQQs at z ^ 4 detected by Chan¬ 
dra with more than 2 counts, we have derived X-ray 
spectral constraints for the quasar population at high 
redshift. A power law with r=1.93±0.09 is a good pa¬ 
rameterization of the quasar spectra in the rest-frame 
Ri 2^0 keV band, similar to what is found for the 
hard X-ray continua of low- and intermediate-redshift 
quasars. Our results are consistent with a lack of signif¬ 
icant spectral evolution of quasar X-ray continua over 
cosmic time. A partial correlation analysis of the pho¬ 
ton index vs. redshift/luminosity for a uniform sample 
of quasars over a broad range of redshift and luminosity 
is required to address this issue further. 

The tight constraint on the absorption intrinsic to the 
quasars. Ah 5.0 x 10^* cm^^, indicates the lack of 
widespread obscuration. 
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TABLE 1 

Chandra OBSERVATIONS OF HIGH-REDSHIFT QUASARS 


Object 

Name 

Z 

Optical 

(CC2000) 

Optical 

(52000) 

Aopt-x** 

(arcsec) 

X-ray 

Obs. Date 

Exp. Time*’ 

(ks) 

Ref. 

BR 0331-1622 

4.36 

03 34 13.4 

-16 12 04.8 

0.4 

2003 Jun 17 

4.67 

1,2 

BR 0353-3820 

4.55 

03 55 04.9 

-38 11 42.3 

0.3 

2003 Sep 10 

4.06 

1,2 

BR 0418-5723 

4.46 

04 19 50.9 

-57 16 13.1 

0.7 

2003 Jun 15 

3.99 

1,2 

BR 0424-2209 

4.32 

04 26 10.3 

-22 02 17.6 

0.4 

2002 Dec 14 

4.67 

1,2 

PSS 0747+4434 

4.43 

07 47 49.7 

+44 34 16.9 

0.5 

2002 Dec 17 

4.54 

1 

PSS 1058+1245 

4.33 

10 58 58.4 

+ 12 45 54.8 

0.1 

2003 Mar 2-3 

5.05 

3 

BR 1117-1329 

3.96 

11 20 10.3 

-13 46 25.7 

0.4 

2003 Jan 28 

4.70 

4 

PSS 1506+5220 

4.18 

15 06 54.5 

+52 20 05.2 


2003 Mar 14 

4.86 

3 

PSS 1646+5514 

4.04 

16 46 56.4 

+55 14 46.4 

0.3 

2003 Sep 10 

4.85 

1 

BR 2213-6729 

4.47 

22 16 51.9 

-67 14 43.4 

0.4 

2003 Mar 10 

4.90 

1,2 

PSS 2344+0342 

4.24 

23 44 03.2 

+03 42 25.5 


2003 Nov 20 

5.10 

1 


References. — (1) Peroux et al. 2001; (2) Storrie-Lombardi et al. 2001; (3) Omont et al. 2001; (4) Storrie-Lombardi et al. 1996. 

Note. — The optical positions of the quasars have been obtained using SEXTRACTOR (Berlin & Amouts 1996) on the DPOSS2 images, while the X-ray 
positions have been obtained with WAVDETECT. 

^Distance between the optical and X-ray positions; the dots indicate no X-ray detection. 

'^The Chandra exposure time has been corrected for detector dead time. 
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TABLE2 

X-RAY Counts, Band Ratios, and Effective Photon Indices 




X-ray Counts^ 




Object 

[0.3-0.5 keV] 

[0.5-2 keV] [2-8 keV] 

[0.5-8 keV] 

Band Ratio’’ 

pb 


BR 0331-1622 

< 3.0 

8 Q+‘^-^ 
o.y_2.9 

o 0+3.2 

:).v_i 9 

i2-+3:5 

0 '^4+0-31 

1 6+°-'' 
’■'’-0.6 

BR 0353-3820 

4.0+5:^ 

45 9+'^-^ 

10.91^:^ 

56.7+*'® 

fk 1 O+0.08 
^•^^-0.06 

2.0+0.3 

BR 0418-5723 

<3.0 

'■^-2.6 

<6.4 

9.01^-’ 

<0.70 

>0.8 

BR 0424-2209 

<4.8 

io.o+;^'5 

< 6.4 

n.9tti 

<0.49 

> 1.2 

PSS 0747+4434 

<4.8 

6.0+ij 

<3.0 

e.otli 

<0.39 

> 1.6 

PSS 1058+1245 

<4.8 


<6.4 

s.otli 

< 1.64 

>0.1 

BR 1117-1329 

<3.0 


<3.0 

l.Otfl 

< 1.15 

>0.5 

PSS 1506+5220 

<3.0 

<4.8 

<4.8 

<6.4 



PSS 1646+5514 

<3.0 

e.otli 

<3.0 

e.otli 

<0.38 

> 1.4 

BR 2213-6729 

<3.0 

n.9tli 

3 

.’■’’-1.6 

20.9tll 

0 13+^1^ 

9 O+0.7 
^•■^-0.6 

PSS 2344+0342 

<3.0 

<3.0 

<3.0 

<3.0 




^ Errors on the X-ray counts were computed according to Tables 1 and 2 of Gehrels (1986) and correspond to the iG level; these were calculated using Poisson 
statistics. The upper limits are at the 95% confidence level and were computed according to Kraft, Burrows, & Nousek (1991). For sake of clarity, upper limits 
of 3.0, 4.8, and 6.4 indicate that 0, 1, and 2 X-ray counts, respectively, have been found within an extraction region of radius 2” centered on the optical position 
of the quasar (considering the background within this source extraction region negligible). 

^ We calculated errors at the 1 g level for the band ratio and effective photon index following the “numerical method” described in § 1.7.3 of Lyons (1991); 
this avoids the failure of the standard approximate variance formula when the number of counts is small (see § 2.4.5 of Eadie et al. 1971). The band ratios and 
photon indices have been obtained applying the correction required to account for the ACIS quantum efficiency decay at low energy. 
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TABLE 3 

Optical, X-ray, and Radio Properties oe the Sample of High-Redshift Quasars 


Object 

(1) 

1Vh“ 

(2) 

^®1450{l+z) 

(3) 

fisoo^ 

(4) 

log (vLv) 2500 
(5) 

Mb 

(6) 

Count rate‘s 
(7) 

f/ 

(8) 

fl kev'^ 

(9) 

log(vLv)2 keV 
(10) 

log(7.2-10 kev)^ 
(11) 

a ^ 

(12) 

r'^ 

(13) 

BR 0331-1622 

5.99 

17.7 

4.64 

47.3 

-29.3 

i.9o;!J:** 

8 

°-^-2.8 

6.81 

45.1 

45.3 

1 q#;+0-09 
^•°°-0.08 

<2.9‘ 

BR 0353-3820 

1.51 

17.9 

3.86 

47.2 

-29.2 


44 7+2.6 

37.04 

45.8 

46.0 

— 1 34+0-07 

5.2‘ 

BR 0418-5723 

1.72 

17.8 

4.24 

47.3 

-29.2 

1 

^•^°-0.66 


5.69 

45.0 

45.2 

1 ©"T+O-O^ 


BR 0424-2209 

2.51 

18.0 

3.52 

47.2 

-29.0 

2 14+0-92 
^•^^-0.66 


6.84 

45.1 

45.3 

-1.81+0.08 

<3.8‘ 

PSS 0747+4434 

5.25 

18.4 

2.44 

47.0 

-28.6 

1 t'->+0.79 
^••^^-0.53 

5 7+3-^ 
■^■'-2.2 

4.66 

44.9 

45.1 

-1.81+0.10 

< 1.5J 

PSS 1058+1245 

2.05 

17.7 

4.64 

47.3 

-29.3 

n co+0.58 
^•^^-0.31 

2Atfj 

1.88 

44.5 

44.7 

— 2 07+0-13 

<0.8i 

BR 1117-1329 

5.05 

18.3 

2.67 

47.0 

-28.5 

0 43+0-57 

LStji 

1.37 

44.3 

44.5 

—2 03“*"^'^^ 

<4.T 

PSS 1506+5220 

1.68 

18.3 

2.67 

47.0 

-28.6 

<0.99 

<3.9 

<3.02 

<44.7 

<44.9 

< -1.90 

< 1.41 

PSS 1646+5514 

2.34 

17.4 

6.12 

47.4 

-29.5 

^•^^-0.50 

^•^-2.0 

3.78 

44.8 

45.0 

-2.00+0.10 

< o.E 

BR 2213-6729 

2.73 

18.5 

2.22 

47.0 

-28.6 

•^•‘^^-0.86 

14.9l^'^ 

12.17 

45.3 

45.5 

-1.64+0.07 


PSS 2344+0342 

5.56 

18.2 

2.93 

47.1 

-28.8 

<0.59 

<2.6 

<2.05 

<44.5 

<44.7 

< -1.98 

<4.5‘ 


Note. — Luminosities are computed using Hq=10 km s ^ Mpc ^ , ^2m=0-3, and 
^From Dickey & Lockman (1990) in units of 10^^ cm“^ . 

*^Rest-frame 2500 A flux density in units of 10“^^ erg cm“^ s“^ Hz“^ 

‘^Observed count rate computed in the 0.5-2 keV band in units of 10“^ counts s“^ 

^Galactic absorption-corrected flux in the observed 0.5-2 keV band in units of 10“^^ erg cm“^ s“^. These fluxes and the following X-ray parameters have been corrected for the ACIS quantum efficiency decay at 
low energy. 

^Rest-frame 2 keV flux density in units of 10~^^ erg cm“^ s“^ Hz“^. 

^Rest-frame 2-10 keV luminosity in units of erg s“^. 

^Errors have been computed following the “numerical method” described in § 1.7.3 of Lyons (1991); both the statistical uncertainties on the X-ray count rates and the effects of the observed ranges of the X-ray and 
optical continuum shapes have been taken into account (see §3 of VO 1 for details). 

*^Radio-loudness parameter, defined as R = fs cnz/f^QQ ^ (rest frame; e.g., Kellermann et al. 1989). The rest-frame 5 GHz flux density is computed from the observed 1.4 GHz flux density (mainly taken from the 
FIRST and the NVSS) assuming a radio power-law slope of a = —0.8, with /v c>c v“. For two quasars, BR 0418—5723 and BR 2213—6729, no radio measurements are available. 

H.4 GHz flux density from FIRST (Becker, White, & Helfand 1995) . 
j 1.4 GHz flux density from NVSS (Condon et al. 1998). 

*^1.4 GHz flux density from Carilli et al. (2001). 
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TABLE 4 

Results from Joint Spectral Fitting 


Quasars used 

Number of 



X-ray 


Wh,z“ 

in the fitting 

Quasars 

Median z 

Median Mb 

Counts 

r 

(cm-2) 

All RQQs’’ 

48 

4.43 

-28.54 

872 

1 na+O-lO 

< 4.95 X 10^* 

All but the 3 highest S/N QSOs 

45 

4.43 

-28.48 

674 

1.92±0.10 

< 7.48 X 10^* 

Most optically luminous half of the sample 

24 

4.30 

-29.06 

563 


< 6.50 X 10^* 

Least optically luminous half of the sample 

24 

4.60 

-28.25 

309 

1.88+'>:1* 

< 1.19 X 10^2 

Highest redshift half of the sample 

24 

4.83 

-28.10 

373 

1 

^•^^-0.16 

< 1.76x 10^2 

Lowest redshift half of the sample 

24 

4.20 

-28.75 

499 

1.92±0.12 

<4.86x 10^‘ 


^Column density in the quasar rest frame. 

*^The slightly different spectral-fit result reported in §6 (F — 1.93 ±0.09) is referred to the power-law model without absorption. 
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Fig. 1.— Full-band (0.5-8 keV) images of the quasars observed by Chandra and presented in this paper. The panels are 98'^ x 98^' for all of the images, except for the source 
detected with two counts (BR 1117—1329; the two counts fall in the same pixel) and those without X-ray detections (PSS 1506-1-5220 and PSS 2344-1-0342); for these three sources, 
the panels are 49'' x 49". North is up, and East to the left. The images of the sources with > 2 counts have been adaptively smoothed at the 2a level; the optical positions of the quasars 
are marked by crosses. For the three sources with < 2 counts, the raw (un-smoothed) images are shown; in these cases, 2"-radius circles around the optical positions of the quasars are 
shown. We note that, although two photons fall within the 2"-radius circle centered on the position of PSS 1506-1-5220, this is not considered a detection (see §2.1). 
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Wavelength (A) 


Fig. 2.— The spectrum of PSS 1506+5220 obtained with the Low Resolution Spectrograph of the Hobby-Eberly Telescope. The data have a spectral resolution 
of 1100. The two prominent broad absorption troughs by Si IV and CIV are labelled, as are the strong Mg II absorption doublet at z=1.4711 and the expected 
location of the C IV emission line at two redshifts (see §3.2 for details). 
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Fig. 3.— Observed-frame, Galactic absorption-corrected 0.5-2 keV flux versus ABi 45 o(i-i- 2 ) magnitude for z >3.96 radio-quiet AGN, mainly quasars. The 
objects presented in this paper are shown as large filled triangles (X-ray detections) and large downward-pointing arrows (upper limits); RQQs from previous 
X-ray observations (Kaspi, Brandt, & Schneider 2000; VOl; Brandt et al. 2002; Bechtold et al. 2003; V03a; V03b) are shown as filled circles (X-ray detections) 
and small downward-pointing arrows (upper limits). AGN detected by X-ray surveys (Schneider et al. 1998; Silverman et al. 2002; V02; Castander et al. 2003; 
Treister et al. 2004) are shown as open squares. The slanted dashed lines show the aox=—1-5 and aox=—1-8 loci at z=4.3 (the average redshift of the present 
sample); the dotted line shows the best-fit correlation reported in §5 for the sample of high-redshift, optically selected RQQs observed by Chandra. 




Number of Quasars Number of Quasars 
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Redshift 



Redshift 




Exposure Time (ks) Redshift 


Fig. 4. — Relevant information for the 48 quasars used in the X-ray spectral analysis: (a) histogram of redshift; (b) redshift vs. absolute B-band magnitude 
(the sizes of the symbols increase with the number of counts); (c) histogram of the exposure times (corrected for detector dead time); (d) redshift vs. number of 
counts and relative errors (computed following Gehrels 1986). 
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CQ 




Observed Energy (keV) 


Fig. 5.— (a) 68, 90, and 99% confidence regions for the photon index vs. intrinsic column density derived from joint X-ray spectral fitting of the 48 RQQs at z >3.99 detected 
by Chandra with more than 2 counts. The average redshift of the sample is z=4.57, while its median redshift is 2=4.43. The solid contours have been obtained using all of the source 
counts, while the dashed contours show the results obtained in the rest-frame 2.2^0 keV band common to all of the quasars of the present sample, (b) Combined spectrum of the 48 
RQQs at z >3.99 fitted with a power law and Galactic absorption (top panel) and data-to-model residuals in units of <5 (bottom panel). 
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Mb 


Fig. 6.— (a) Hard X-ray photon index vs. redshift for a compilation of RQQs. The circles at z < 1.8 indicate RQQs from the Bright Quasar Survey (PG) 
analyzed by Piconcelli et al. (2005), while the diamonds indicate the z 1.8 — 2.5 RQQs from Vignali et al. (2001b) sample. At higher redshift, the data points 
indicate joint spectral htting (triangles) and single-object spectroscopy (squares); numbers in the data points correspond to the citations given in the figure legend, 
(b) Hard X-ray photon index vs. Mb for the same quasars. 







Fig. 7.— ttox distribution for the optically selected RQQs observed by Chandra at z > 3.96; note that here we have included also the two-count source 
BR 1117—1329 at z=3.96 which was not used in the joint X-ray spectral fitting. The shaded area indicates the BALQSOs; the areas with arrows indicate the 
95% confidence upper limits. The long-dashed (dotted) line indicates the average ttox (using ASURV) for high-redshift non-BAL RQQs (BALQSOs) observed by 
Chandra. The typical uncertainty in the ttox measurements is ^ 0.15. 
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TABLE A5 

Properties oe z > 4 SDSS DR2 Quasars in ROSAT Fields 


Object 

Z 

Afil450(l+z) 

Mb 



Ot-ox 

R'’ 

SDSS 004054.65-091526.79 

4.973 

19.3 

-27.6 

< 

3.24 

< 

-1.34 

<5.6 

SDSS 100645.60+462717.25 

4.440 

20.1 

-26.7 

< 

3.02 

< 

-1.29 

< 136.7 

SDSS 101053.00+531144.82 

4.506 

20.1 

-26.7 

< 

1.14 

< 

-1.45 

< 10.1 

SDSS 123735.47+642936.00 

4.334 

20.1 

-26.6 

< 

4.22 

< 

-1.23 

<32.7 

SDSS 140146.53+024434.72 

4.375 

18.7 

-28.0 

< 

2.36 

< 

-1.54 

<2.8 

SDSS 150730.63+553710.83 

4.499 

20.2 

-26.6 

< 

1.71 

< 

-1.36 

< 11.1 

SDSS 213243.26+010633.91 

4.032 

20.2 

-26.4 

< 

3.65 

< 

-1.25 

< 10.5 


“Galactic absorption-corrected flux in the observed 0.5-2 keV band in units of 10 erg cm ^ s ’ 
'’Radio-loudness parameter. 


APPENDIX 

NEW ROSAT UPPER LIMITS FOR Z >4 SDSS DR2 QUASARS 

Here we report some additional X-ray flux upper limits for the Sloan Digital Sky Survey Data Release 2 (SDSS DR2; Abazajian 
et al. 2004) quasars at z > 4. These were obtained with ROSAT PSPC and HRI using only the inner w 20' of these instruments. 
The 3o upper limits (see Table Al) have been derived in a similar way to those shown in Table 4 of VOl using the SOSTA 
command in the XIMAGE package (Giommi et al. 1992) Version 4.2. The AB 145011 + 7 ! magnitudes have been derived using the 
method described in §2.2 of VBS03. 
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TABLEB6 

Optical, X-ray, and Radio Properties oe z > 4 Quasars in Archival Chandra Observations 


Object 

Z 

^■Sl450(I+z) 

f 

''2500 A 

\os(L 0^ 

2500 a' 

Mb 

/o.5-2 keV 

h keV 

log (7-2 kev)'’ 

log(72-10kev) 

Ctox 

R 

PSS 0134+3307 

4.53 

18.8 

1.69 

31.80 

-28.3 

11.5 

9.45 

27.55 

45.44 

-1.63 

< 8.1 

PSS 0808+5215 

4.44 

18.2 

2.93 

32.03 

-28.8 

3.5 

2.81 

27.01 

44.90 

-1.93 

0.2 

BR 1202-0725 

4.70 

18.0 

3.52 

32.14 

-29.1 

8.5 

7.21 

27.45 

45.34 

-1.80 

< 1.3 

BR 1600+0724 

4.38 

19.3 

1.06 

31.58 

-21.1 

6.8 

5.43 

27.28 

45.21 

-1.65 

<3.4 

BR 2235-030E 

4.25 

18.9 

1.54 

31.72 

-28.1 

<2.6 

<2.03 

< 26.84 

< 44.73 

< -1.87 

< 8.5 

PSS 2322+1944 

4.12 

17.9 

3.86 

32.09 

-29.0 

10.4 

7.92 

27.41 

45.30 

-1.80 

0.2 


Note. — Units are the same as in Table 3. 

^Rest-frame 2500 A luminosity density (erg s~^ ). 

Rest-frame 2 keV luminosity density (erg s~ ^ Hz” ^). 
‘^BALQSO. 


OPTICAL, X-RAY, AND RADIO PROPERTIES OE Z > 4 QUASARS IN ARCHIVAL Chandra OBSERVATIONS 

Here we report the main optical. X-ray, and radio properties for z > 4 quasars with unpublished Chandra observations. Note 
that BR 1202—0725 was detected also by ROSAT (Kaspi, Brandt, & Schneider 2000). 








